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ABSTRACT 

A random, hydrogen-free, assembly of microscopic sp^ carbon chips, forming a 
macroscopically homogeneous and isotropic sohd, is proposed as a model car- 
rier for the UV interstellar extinction band . The validity of this model is based 
on the calculation of the Bruggeman average dielectric function of a mixture 
of the known parallel and perpendicular dielectric functions of graphite. The 
TT absorption feature of Rayleigh-sized spheres of this mixture falls near 4.6 
//m^^ (2175 A ), but its width is 1.5 /im^^, somewhat larger than the astro- 
nomically observed average, 1 fim^^. This is confirmed by measurements of 
the reflectance of an industrial material, polycrystalline graphite. A better fit 
to the IS feature position and width is obtained with a hypothetical mate- 
rial, having the same dielectric functions as natural graphite, except for less 
extended wings of the n resonance. Physically, this could result from changes 
in the electronic band structure due to previous thermal histories. On this 
model, the Frolich feature central wavelength depends only on the tt reso- 
nance frequency, while its width depends only on the damping constant of 
the same resonance. This explains the range of observed feature widths at 
constant feature wavelength. 

Key words: Stars:Astrochcmistry; interstellar matter; ISM: dust; molecular 
processes. 
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1 INTRODUCTION 

The numerous carbon models that have been proposed to mimic the 2175 A (or 4.6 /im~^) 
interstellar (IS) extinction band can roughly be divided into two generic groups respectively 
characterized by their ordered or disordered solid structure. Typical of the first group is 
the distribution of randomly oriented pure graphite spheres (see Draine and Lee (119841) ). 
FuUerenes and bucky onions (BO) were later proposed (e.g. De Heer and Ugarte (119931) . 
Wada and Tokunaga fl2006p ). These recently returned to the forefront, with the laboratory 
isolation of single BOs by Chhowalla et al. (I2003P : their spectral extinction was measured 
and shown to exhibit a band at 4.55±0.1/im~^ whose profile, when adequately smoothed, 
fits the "Drude-like" profile proposed by Fitzpatrick and Massa (119861) . This experimental 
result is backed by theoretical analyses, one of the latter being that of Ruiz et al. (120051) . 
who demonstrated a very good agreement of their result with the measurement. 

The main difficulties faced by such "ordered" models are (see Draine and Lee (119841) . 
Draine and Malhotra (fT993l) ): 

a) the absence of a known process by which the FWHM of the band could vary by a 
factor of ~2, as observed, even as its peak and shape remain unaltered; 

b) in particular, the measured range of FWHMs of BOs is 1.2-1.6 fim~^, outside the 
observed range (0.7 to 1.3 /xm"^); it is not clear how BOs could be altered to yield FWHMs 
in the required range; 

c) the appearence, in the spectrum of the model, of several features that are not observed 
in the sky, e.g. the four IR features of the Cgo molecule (see Huffman (119910 ): 

d) the absence of efficient IS processes capable of producing enough fuUerenes or BOs, 
which are known to be hard to obtain, even in specialized laboratories (see Chhowalla et al. 
( I2003P ). The same applies to pure graphite particles. 

These difficulties prompted research in the opposite direction, namely amorphous carbon. 
This ill-defined term covers a large family of materials, hypothetized or indeed known in the 
laboratory (see Robertson (I1986P ). and characterized by an array of vr and a bonds randomly 
distributed in space, on the microscopic scale. This is not to be confused with random 
macroscopic clusters of graphitic particles (e.g. Draine ( I1988P ). or composite grains (e.g. 
Mathis and Whiffen (I1989P ). One of the first disordered carbon model for the 2175 A band 
was proposed by Sakata et al. (I1983p ) who synthesized a quenched carbonaceous composite 
material (QCC), displaying a UV absorption near 2000 A . A number of workers followed 
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suit, (for references, see Bussoletti et al. ( I1995p . Common to these workers is, apparently, 
the idea exphcited by Hecht fll986p . based on the work of Donn (119680 and Nuth (11985P 
on kerogen-hke models of IS grains, and the evolutionary trend of such earthly materials 
upon heating, as observed, in particular, by Fink et al. (I1984p : it is suggested that, in this 
graphitization process, the H atoms they initially contain are progressively expelled while 
the essentially C residual atoms are spatially redistributed, allowing somehow the 2175 A 
feature to emerge. 

Pursuing this effort, we return here to our polycystalline graphite (PG) model (Papoular 
et al. (I1993p . using a new approach and bringing new evidence in favour of such a model. 
In Sec. 2, we first show on general grounds, how the position and width of the band can be 
disconnected for an amorphous material, and what their governing factors are. 

As an embodiment of this notion, we envision in Sec. 3 an isotropic homogeneous mate- 
rial made up of randomly oriented, atomic-scale, graphitic lumps; its dielectric function is 
computed, according to the Bruggeman rule, by mixing the lumps in adequate concentra- 
tions. This is shown, in the Rayleigh approximation of small nearly spherical grains, to yield 
a strong band at 20 A away from the observed IS wavelength, with a FWHM of 1.6 /im"^. 

We then tackle the problem of reducing the band width to the IS standards without 
letting the central wavelength drift excessively: this is done by varying the dielectric prop- 
erties of the constituent "graphite" only to a very limited extent. In this process, the band 
retains a shape that reasonably fits the "Drude-like" profile and its intensity satisfies the IS 
C- abundance constraint. 

Section 4 compares these results with our previous measurements of man-made PG (Pa- 
poular et al. (119931) ): it is shown that the latter are adequately mimicked by our model 
associated with common graphite, lending more credibility to the "amorphous" model. This 
provides a link to models using less graphitized constituent materials, which may help ex- 
plaining the many observed "red" outlier bands. Finally, Section 5 discusses in more detail 
the compatibility of the model with observational constraints. 

2 A PROPERTY OF THE FROHLICH RESONANCE 

The weak scattering apparently included in the UV band extinction, elicits the assumption 
of small particles for the band carrier, and, hence, the application of the Rayleigh approx- 
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imation. As a corollary, surface effects cannot be neglected, which entails the use of the 
Frohlich expression for the extinction efficiency of a small sphere of radius a: 
r., 247r 

where ei,2 are, respectively, the real and imaginary parts of the dielectric function of 
the assumed homogeneous material; L is the shape parameter of the supposedly ellipsoidal 
grain, 1/3 in the case of a sphere. We assume throughout that scattering is negligible so that 
there is no need to distinguish between absorption and extinction efficiencies, both being 
designated by Q. 

In order to quantitatively make our point, we focus on a spectral region where the 
dielectric function of the material can be represented by a Lorentzian oscillator (Bohren and 
Huffman (119831) ) with resonant frequency cuo; plasma frequency ojp and damping constant 7 
(assumed small): 

= 1 + - + h/2)- = (.„ - + (V2)- 

Plots of ei^2 and Q/a for different values of the Lorentzian parameters show that, there 
being spectral dispersion, the maximum of Q does not occur when ei = 1 — L~^, which is 
true only if 7/^79 is vanishingly small; ei need not even fall below 1 — L^^ for a Frohlich 
resonance to occur. In the case of graphite, contrary to crystals like SiC, for instance, the 
excursion from this equality is considerable and depends on the slope of e vs frequency. 
Inspection of the plots suggests instead that 

ei + 2 = €2, (3) 

The numerical fact that this relation is approximately verified near the maximum of 
Q/a can be understood by noting that, in the vicinity of our resonance, ei^2 both vary, but 
in opposite senses. Since they appear in a sum of two terms in the denominator of Q/a (eq. 
1) it seems reasonable to assume, as an approximation, that the maximum (at Um) occurs 
when the two terms are equal. 

Relation 3 is the more closely verified the smaller 7 is with respect to Uq and Up. If it is 
not too large, one readily concludes that 

UJm^ UJ0 + 77^. (4) 

Quo 

For an ellipsoid in general, the factor 6 in the denominator must be replaced with 2/L, 
meaning that the small variations of the Frolich resonance frequency are nearly linear in L. 
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Note the absence of 7 in this expression. Now to evaluate the FWHM of the resonance. By 
definition, when uJi/2 = ciJm±FWHM/2, 

^/ ^ / 247r 1 , , 

Q/a = Qmaxl'2^a = ^• 

Rewriting this expression in terms of the 3 Lorentzian parameters, and assuming again 
that 7 << tuo.p! one obtains 

cui/2 ~ cu^ ± ^ ; FWHM = ^87. (6) 

Thus, the sharpness of the Frohhch resonance increases hnearly with the steepness of the 
Lorentzian (blue) wing. Here, only 7 appears. That the dependencies of FWHM and peak 
frequency are distinct is not specific of the Lorentzian but results from the Kramers-Kronig 
link between the real and imaginary parts of the dielectric function. This result will help, 
below, tailoring the width of the model UV extinction band. 



3 THE MODEL CARRIER 

It is now widely recognized that IS carbon grains cannot be pure graphene chunks, but also 
include "impurities" (mainly hydrogen, oxygen and nitrogen) which destroy the beautiful 
symmetry of the structure, thus determining the direction of the subsequent evolutionary 
trend. The energy absorbed by the grain in the harsh IS medium progressively removes 
the volatile impurities, which are less tightly bound than carbon. Starting from a disordered 
state, this process can only add disorder to the structure. Our model carrier is in the ultimate 
stage of this evolution, where, in order to obtain a vr-type resonance, we have to assume that 
all adatoms have been expelled, and the carbon bonds are mainly of the sp^ type, like in 
polycrystalline graphite or glassy carbon (Taft and Philipp ( I1965p . Robertson (I1986P ). More 
precisely, this can be construed as an aggregate of small stacks of graphene planes which, 
instead of piling up one upon the other, normally to their common direction, would lie at a 
random angle to one another. This arrangement is inspired by the structure first described 
by Franklin (I195ip for "non-graphitizing" , later called glassy carbon. 

The size of each stack is too small for the stacks to be considered as ellipsoids to which 
Frohlich's theory may be applied. Instead, we consider that the constituent material of 
the grain is microscopically homogeneous, with electromagnetic properties resulting from 
an adequate mix of two different materials. Indeed, along any given direction, each stack is 
viewed either approximately parallel or approximately perpendicular to the graphene planes. 
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It is assumed here that the situation is akin to that of a mixture of two homogeneous and 
isotropic materials, whose dielectric functions are, respectively, e || and e± of graphite, 
with corresponding concentrations of 1/3 and 2/3. It is usually considered that, for such 
mixtures, the Bruggeman mixing formula is a good choice among the host of proposed 
formulae (Bohren and Huffman (I1983p ). In the present instance, we write this formula as 

<^ II +^ee// e± + leeff 

where f\\ = 1/3 and ee// is the required isotropic, spatially uniform dielectric function. 
The factor 2 in the denominators is the particular value of a function which depends on 
the geometrical factor, L, defining an ellipsoid, and is valid for spheres {L = 1/3). This 
procedure has recently been rediscussed (Lucarini et al. (I20U5P ). and shown to satisfy the 
Kramers-Kronig causality relationship. 

It must be stressed that, for complex dielectric functions, the Bruggeman equation admits 
of several solutions, among which a physically sensible choice must be made. A copy of the 
code we have written can be obtained upon request. For e ||, we used the tabulated values in 
Draine (I1985P : for e±, we used both the original values of Taft and Philipp (I1965P and those 
adjusted by Draine and Lee fll984p to deliver the right UV band according to their model 
of randomly oriented pure graphite spheres. Figure 1 shows that the Draine and Lee values 
lead to a somewhat better overall band shape, when compared with the fittest "Drude" 
profile of Fitzpatrick and Massa (I1986p . 

Both features peak at 2200 A , near the observed wavelength; however they are both 
wider (1.5 yum^^) than the astronomically observed average, 1 fim^^. Some tailoring of the 
dielectric properties is therefore in order, and Sec. 2 may help in doing this. 

According to eq. 4, a reduction of only ~ 3% in the vr plasma frequency is enough to 
bring the extinction peak frequency to the desired value. Such a change is conceivable in 
space, as Up is intimately linked to the electrical conductivity, which is very sensitive to the 
heat treatment previously suffered by the material (Robertson (I1986p ). Alternatively, the 
feature location may be controlled through e ||, which is notoriously subject to experimental 
errors of measurement (see Jellison et al. (I2007P ) and was recently shown to strongly depend 
on the distance between graphene planes (Fei et al. (120081) . Marinopoulos et al. (120021) ). 
Thus, calculations show that changing its value by a factor of 1/2 or 2 shifts the feature to 
2070 and 2210 A , respectively; the required adjustment may thus fall within the limits of 
measurement errors and structural variations. 
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Figure 1. The UV extinction feature deduced from the Bruggeman mixing formula applied to e || and e± (respectively 1/3, 
2/3). Squares: t \\ and e± from Draine H1985|l : pluses: e± from Taft and Philipp H1965|l : line: "Drude" profile with Ao=2200 A 
, 7 = 1.5 iira~^ . 



As for the band width, eq. 6 shows that it can be reduced by steepening the blue wing of 
the vr resonance. Now, sohd state physicists have shown that the Lorentzian representation 
of the dielectric function is suited to cases where a single decay channel for excitation exists 
or is dominant. 

On the other hand, the Gaussian representation, which results in a quicker evanescence 
of the wings, is known to be preferable when the decay is primarily associated with a large 
number of weak scattering centers (Garland et al. (11988P and references therein; Franke et 
al. (119961) ). which is likely the case for our disordered model. It was used, for instance, by 
Brendel and Bormann (I1992p for amorphous solids. 

One combination of Lorentzian and Gaussian representations is the so-called pseudo- 
Gaussian (Kim (I1992p ). obtained by replacing the constant 7 in the Lorentzian by a function 
of the form 




7o 



(8) 
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Table 1. Oscillator parameters 



Polarization 




e 






€± 




Resonance 






7 




UJp 


7 


a 


1.7 10^^ 


1.8 10^^ 


SIO^^ 


2.1710^^ 


2.8 10^® 


410^5 


TV 


5.9410^^ 


8.5 10^^ 


410^5 


6.731015 


1.12 lOl** 


77r± 



All values in rad.s ■yTr± from eq. 8. 



This was found to be helpful, in particular, for terrestrial graphite (Djurisic and Li 
(I1999P ). Here, we used the latter device to represent the 7r± resonance, and Lorentzians for 
the other 3 resonances, e// and e_L are each the sum of the corresponding a and vr functions. 
Table 1 lists all the constants used to compute 3 band spectra, covering a representative 
range of FWHM widths; the last entry, 7,^^, is given by eq. 8, where 70 = 2 lO"*^^ rad.s^^ and 
a was given 3 values, successively: 1/4, 1, 4; only a was changed in the process, not uq, nor 
Up, so as to limit the excursion from natural graphite properties to a minimum. Figure 2 
illustrates the results. 

In Fig. 2, the upper and lower graphs, respectively, plot ei and 62 for the 3 values of 
a: 1/4, 1, 4, as well as the corresponding curves for a Bruggeman mixture e\\/3 + 2e^/3, 
taken from Draine fll985p . Note the limited extent of differences between the synthesized 
dielectric functions of the 3 cases, in the UV band window, 3-6 fim~^. Below 4 fim~^, 
most of this difference is due to the simplistic representation of the material by only 2 
pairs of Lorentzian oscillators. The middle graph represents the corresponding extinction 
efficiencies in the Rayleigh limit. Fitting a "Drude" function to each delivers the width and 
corresponding central wavelength 

7=1.6, 1.35, 1 fim~^ 

Ao=2150, 2150, 2190 A . 

Figure 3 illustrates the trends of 7 and Xmax as a varies. In agreement with Fig. 2 of 
Fitzpatrick and Massa (1986), Xmax remains within 1.1 % of 2175 A . 
Figure 4 shows the Drude fit for the narrowest feature. 

Even though this fit may be considered as acceptable, it is limited to a small wavenumber 
interval. Based on the astronomical observations in this interval alone, one cannot infer 
the overall properties of the carrier. However, the variations of the dielectric properties, 
required to cover the diversity of the observed features, are limited, and the difference is 
small between these properties and those of "astronomical" graphite. More quantitatively, 
the effective electron number as a function of frequency, as deduced from the Sum Rule (see 
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Figure 2. Middle graph: Q/a (A of a sphere of radius a, made of the model carrier (Sec. 3), in the spectral window 
of interest, ei and £2 in eq. 1 were obtained from the Bruggeman averaging formula (7), successively applied to 3 pairs of 
Lorentzian-like e ||, J- (see text and Table 1). Only one parameter, a, characterizing the width of the ttI. resonance varies from 
case to case; a=0.25 (line), 1 (dashes), 4 (dotted), respectively. The extinction feature grows and narrows as a increases. Upper 
and lower graphs: the 3 corresponding ei and 62 ■ For comparison, the heavy dots represent the average dielectric functions 
obtained from the Bruggeman formula by inserting the dielectric functions tabulated by Draine II1985I I for e || and el. of 
"astronomical" graphite. 
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Figure 3. The trends of 7 and Aq as a varies (eq. 4 and 6). In agreement with Fig. 2 of Fitzpatrick and Massa (1986), A„ 
remains within 1.1 % of 2175 A . 



Taft and Philipp (119651) ) for the cases considered in Fig. 2 is hardly distinguishable from 
that of astronomical graphite. A reasonable surmise is, therefore, that the latter may be 
taken as a first approximation to the functions of our model carrier. 



4 A LABORATORY MATERIAL EMBODIMENT: POLYCRYSTALLINE 
GRAPHITE (PG) 

High purity PG is produced industrially for various purposes. It is not to be confused with 
HOPG, or Highly Oriented Pyrolytic Graphite, which is also produced under high temper- 
atures and high pressures. HOPG is very close to the perfect natural graphite, while PG 
is perhaps the best yet embodiment of the opposite extreme: a random assembly of mi- 
croscopic sp^ carbon chips pressed together, (also termed "turbostratic crystallites" in the 
Carbon literature; see Mrozowski (119711) ). which we propose for a model band carrier. We 
have measured the normal-incidence reflectance of the polished surface of a pressed pellet 
made from high-purity PG powder (Papoular et al. (119931) ). This provides us now with the 
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Figure 4. "Drude" fit (line) for the narrowest feature of Fig.2 (dots): Ao = 2190 A , A'^ = 4.566 //m-\ 7 = Ifim-^ 

opportunity to test the validity of the theoretical mixing procedure presented above. Figure 
5 displays the measured reflectance of PG, which we now set out to fit with the model of 
Sec. 3. For this purpose, we follow the same procedure, and define 4 oscillators as above. 
Since the material is supposedly pure graphite, we try first the same oscillators as in Table 
1, except that, for a start, we adopt a constant 7 for all 4 oscillators. Trying different values 
of 7(7r _L), all other parameters remaining unchanged, we compute the dielectric function of 
the mixture "e || /3, 2e±/3", using the Bruggeman formula, then derive the corresponding 

optical indices, n and k, using the following formulae: 

n=(M±fI^)i/2 

k = (M±fipZ£i)l/2. 

The normal reflectance is then given by 

p _ (n-l)^+fc^ 
~ (n+l)2+fc2 • 

It is found that the closest fit is obtained with 7 = constant = 2 lO"*^^ rad.s^"*^; but still a 
better fit results from a minor change in Table 1, namely c<jp(7r±) = 7.48 10^^ , which yields 
the line in Fig. 5. The accuracy of the measurement does not warrant a better fit with a 
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Figure 5. Dots: normal incidence reflectance of polished pressed pellet of industrial polycrystalline graphite powder (PG) 
measured at the LURE synchrotron; line: PG model of Sec. 3, with the same parameters as for the narrow feature of Fig. 4, 
except 7(7r±) = 210^^ rad.s"^, constant; uip = 7.4810^^ rad.s"^, giving an approximate fit, which shows the ability of the 
model to mimic a laboratory material approaching the ideal random sp2 structure. 



model including only 2 pairs of oscillators. This one already demonstrates the respective 
roles of the model parameters. Now, using the parameters of Tab. 2, we can also deduce 
the corresponding band intensity, Qabs/O'- This is displayed in Fig. 6, together with the 
accompanying Drude fit: Aq — 2210 A ; 7 = 1.3 //m~^. Thus, this laboratory material 
approaches qualification as a qualitative model, although comparison with fig. 4 shows that 
its feature is much weaker than that of the hypothetical model of Sec. 3. 

A better laboratory analog might be obtained by trying to mimic more closely the evo- 
lutionary trend of carbon grains in the sky. The starting material could be kerogen, whose 
disordered structure certainly makes it a non-graphitizing carbon. A subsequent slow and 
protracted high temperature treatment at low ambient pressure would expel all volatile 
atoms, leaving a nearly pure carbon substance approaching the ideal turbostratic graphite. 

Still another analogy can be found in the insoluble organic matter (lOM) found in car- 
bonaceous chondrites extracted from meteorites like Orgueil and Murchison. A large fraction 
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Figure 6. The UV extinction feature of the laboratory PG of Fig. 5 (dots), together with a "Drude" fit, with Aq = 2210 A , 
7 = 1.3 /xm^^ . 

of the carbon in this material is found in weakly organized aromatic moieties containing 
about 30 to 150 C atoms, whose average structural parameters have recently been measured 
by Transmission Electron Microscopy (see Derenne et al. (120051) ). The average "diameter" 
of these basic structural units corresponds to 2 to 3 aromatic rings per layer, with 2 to 3 
stacked layers per unit. This is reminiscent of the graphitic chunks constituting our model, 
Sec. 3. Interestingly, the average interlayer spacing in measured meteoritic stacks is ~0.48 
nm, much larger than in terrestrial graphite, suggesting differences in dielectric properties. 



5 DISCUSSION 

The properties of the Frohlich resonance (Sec. 2) and the sample calculations of Sec. 3 show 
that the proposed model is able to mimic the UV extinction features, over the observed 
range of position and width, using dielectric properties that differ only slightly from those 
of graphite. The differences occur in a limited spectral range, in the blue wing of the vr res- 
onance, where the measurements of e are themselves affected by considerable experimental 
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uncertainty. Provided the particle size is no larger than ~100 A for the Rayleigh approx- 
imation to be valid, the central wavelength remains within the observational limits, as the 
observed range of width of the astronomical feature is covered by adjusting only the plasma 
frequency and the damping constant of only the 7r± resonance (in fact, only the wings of this 
resonance). We thus concur with the early conclusion of Draine and Malhotra (119931) that 
the causes of variation of the IS feature must lie in the dielectric properties of the carrier. 
These are known to be intimately linked to the details of the solid structure: plasmon-phonon 
interaction, irregularity of the microscopic "graphite bricks" (diameter, interlayer distance, 
etc; see end of Sec. 4), presence of impurities (such as hydrogen), varying strain at interface 
between bricks, size of the latter, etc (see, for instance, Hecht ( I1986p . Sorrell (119901) ). All 
of these depend heavily on the previous thermal history of the material (Franklin (119511) . 
Mrozowski (119711) . Robertson (I1986p ). The detailed investigation of the effects of these fac- 
tors is outside the scope of this paper. However, the particular oscillator parameters listed 
above, together with the laboratory availability of a close analog, PG, should help defining 
further the nature of the proposed band carrier. In the following, we discuss the adequacy 
of the model in more detail. 

5.1 Carbon availability 

Let N{C) and N{H) be the column densities of C and H, respectively, towards a sample 
star, and p, the grain specific gravity. Take 

= 3.1 and Nh = 5 10^-^ Eb-v atoms. cm~^.mag~\ 
and assume that, on average, the optical thickness through the line of sight is t{UV) = 
Ay, then 

N(C) ^ 07 p 
N{H) ^' Q/a- 

For graphite, p ~ 2 g.cm~^, and from Fig. 2, Q/a is of order 10^, assuming small grains 
(less than about 100 A in size), the fraction of cosmic carbon required to be in the model 
grains turns out to be about 54 ppm, or nearly 1/4-1/6 of the available carbon (e.g. Snow 
and Witt (119951) ). Of course the largest fraction must be reserved to less evolved grains. 

5.2 Spurious features 

In agreement with astronomical observations (see Fitzpatrick (120041) ) and due to the absence 
of any alien atom in the ideal model, and to its disordered structure, there should be no 
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features in its spectrum other than the vr and a resonances overlying a continuum; that 
is in fact the case for the measured spectrum of polycrystalhne graphite. It is apparent, 
from the figures above, that the continuum is neghgible in intensity as compared with the 
feature itself. It should therefore make a negligible contribution to that which is observed in 
astronomical spectra. 

On the other hand, the a resonance is much stronger than the vr resonance, and should, 
therefore, be an important component of the far UV rise beyond 6 ^mr^ . Indeed, even when 
the general ascending slope of celestial UV spectra is weak, the local upward bend around 
8 yum"^ remains strong (see the large number of diverse spectra collected in Fitzpatrick 
(I2004p ) . Although available astronomical data do not extend far enough in the UV, they 
carry no sign of incompatibility with the FUV rise in the measured PG spectrum. 

The compact structure of our model drastically reduces the number of surface dangling 
carbon bonds, and therefore, makes it very difficult for adatoms to stick to the grain. This, 
however, remains a subject of further investigation into the possible effects of such adatoms 
on the profile of the vr and a resonances. 

Finally, the perfect isotropy of our model ensures that it causes no polarization of the 
scattered light, provided the grains are spherical or, in any case, poorly aligned. 

5.3 Relation to fullerene models 

In the light of our work, the good fit provided by fullerenes and BOs to the (widest) IS 
features (see Introduction) stems from the fact that, in these particles, the electric dipoles 
associated with the carbon rings are uniformly distributed over the spherical surface, and 
their directions are thus uniformly distributed in angle, exactly as assumed in an ideal 
PG, but on a larger scale. The main problems of the BO model are the electronic spectral 
structure superposed upon the feature, and the large feature width (1.2-1.6 ^mr^\ Chhowalla 
et al. (120031) . Ruiz et al. (l2005|l l 

5.4 Why Bruggeman ? 

The proposed model being a mixture of components, the choice of a mixing formula is 
inescapable. The number of available mixing formulae (see Bohren and Huffman (119831) . 
Stroud (I1975P ) is witness to the theoretical difficulty of the choice. Of course, the arithmetic 
weighted mean is simplest. However, in matters dielectric, it is not clear what is the elec- 
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trie property to be averaged: e.g. resistivities cannot be used in a weighted sum. Moreover, 
common knowledge seems to restrain the use of this ansatz to the inclusion of an impurity 
in an otherwise homogeneous medium. For comparison, we have also used this averaging: 
essentially, it shifts the feature of interest to the red by about 100 A . We have also used the 
Maxwell Garnett mixing formula, which stems from the same basic treatment as Brugge- 
man's but results from a different approximation (Stroud (119750 ). The band obtained in this 
way is shifted, widened and distorted. This does not condemn the mixing formula, but calls 
attention to possible microresonances already observed experimentally (see Barker (119731) ). 
The Bruggeman formula gives a more acceptable feature shape, at the price of a heav- 
ier computation code. Our conclusions regarding the factors governing the feature profile 
and location are not heavily dependent on this choice. It should be interesting to compare 
the results with those of the Discrete Dipole Approximation (DDA) applied to a random 
distribution of dipoles (Draine (11988^ ). 

Another inconvenience with non-linear formulae like Bruggeman's, is the uncertainty as 
to whether they satisfy the Kramers-Kronig relation. We have checked this is the case for a 
mixture of the parallel and perpendicular dielectric functions of graphite (Fig. 1), and, more 
generally, for a sum of two Lorentzian resonators (Sec. 3 and 4). 

5.5 Effects of grain shape 

The Frohlich formula for Qabs and the Bruggeman formula for the average dielectric function 
(see Zeng et al. (119881) ) both include the ellipsoidal characteristic parameter, L, which was 
set above at 1/3, assuming a spherical shape for both the grain and its component subgrains. 
As noted by Bohren and Huffman (119831) . shape effects are no weaker in the Rayleigh ap- 
proximation than for large particles, and the band shape of discs or needles are very different 
from the "Drude" profile. Fortunately, the physics of carbon grain formation and evolution 
apparently does not allow the evolution towards such extreme shapes. Besides, the grain 
surface is certainly very irregular and may not be subjected to the same theoretical treat- 
ment as perfectly smooth ellipsoids. Finally, even though the term "surface mode" is usually 
applied to the Frohlich resonance, it is clear that the whole grain volume is involved, so the 
surface shape may not be so important, after all. In order to substantiate these conjectures, 
we also computed Qabs for values of L symmetrically bracketing 1/3 , from 0.26 to 0.4. It 
was found that, although the peak wavelength of the feature decreases from 2315 to 2140 A 
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, the average feature still very nearly coincides with the curve of Fig. 4. This is due to the 
symmetry of the variations, a consequence of the implicit assumption that the fluctuations 
of L around 1/3 are limited and symmetric. 

6 CONCLUSION 

Small spheroidal particles made up of a homogeneous mixture of randomly oriented sub- 
micron chunks of pure graphite have been shown, above, to display a strong Prolich resonance 
peaking near 2175 A , as does the UV IS extinction feature, the only difference being its 
slightly larger width. A mild, coherent, tailoring of the dielectric function of graphite around 
2000 A suffices to reduce this width to the astronomically observed average. The physical 
justification of such tailoring lies in subtle changes in the damping constant of the 7r± 
resonance, that are likely to result from small alterations of the corresponding electronic 
band structure of graphite in space. These may be due to different grain formation histories 
and/or different IS irradiations and other heating effects. Such variations allow the model 
feature width to cover the range observed in the sky without notably impacting the peak 
wavelength. The Discussion shows that none of the usual observational constraints disagrees 
with this model. 
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Note added in proof: 

The property of separate dependencies of the central wavelength and the FWHM of the 
Frohlich resonance (Sec. 2, eq. 4 and 6) was deduced from the crude approximation in eq. 3. 
If, instead, one goes through the tedious mathematics of differentiating Q{u!) to determine 
the frequencies at the peak and at half-maximum, one obtains the following, more accurate 
relations, to first order in 7/0;, in the relevant range of the parameters: 



UJm^i^O+ ^ (9) 

FWHM = 7 (10) 

3cJo7 



e2m 3 



(11) 



